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Abstract. 

High energy collisions of cosmic ray (CR) nuclei with in- 
terstellar gas are believed to be the mechanism producing 
the majority of CR antiprotons. The distinguishing spectral 
shape with a maximum at 2 GeV and a sharp decrease to- 
wards lower energies makes antiprotons a unique probe of 
the models of particle propagation in the Galaxy and mod- 
ulation in the heliosphere. Besides, accurate calculation of 
the secondary antiproton flux provides a "background" for 
searches for exotic signals from the annihilation of super- 
symmetric particles and primordial black hole evaporation. 
Recently new data with large statistics on the antiproton flux 
have become available which allow for such tests to be per- 
formed. We use our 3D Galactic cosmic ray propagation 
code GALPROP to calculate interstellar propagation in sev- 
eral models. For our best model we make predictions of pro- 
ton and antiproton fluxes near the Earth for different modu- 
lation levels and polarity using a steady-state drift model for 
heliospheric modulation. 



1 Introduction 

Most of the CR antiprotons observed near the Earth are sec- 
ondaries produced in collisions of energetic CR particles with 
interstellar gas (Mitchell et al., 1996). Due to the kinemat- 
ics of the process, the spectrum of antiprotons has a unique 
shape distinguishing it from other cosmic -ray species: it is 
expected to peak at about 2 GeV decreasing sharply towards 
lower energies. In addition to secondary antiprotons there 
are possible sources of primary antiprotons such as candidate 
dark matter particles and evaporating black holes. 

Despite numerous efforts and overall agreement on the 
secondary nature of the majority of CR antiprotons, pub- 
lished estimates of the expected flux significantly differ (see 
e.g. Fig. 3 in Orito et al. 2000). The major sources of uncer- 
tainties are: (i) incomplete knowledge of cross sections for 
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antiproton production, annihilation, and scattering, (ii) pa- 
rameters and models of particle propagation in the Galaxy, 
and (iii) modulation in the heliosphere. While the interstellar 
antiproton flux is affected only by uncertainties in the cross 
sections and propagation models, the final comparison with 
experiment can only be made after correcting for the solar 
modulation. Besides, the spectra of CR nucleons have been 
directly measured only inside the heliosphere while we need 
to know the spectrum outside in interstellar space to compute 
the antiproton production rate correctly. 

We have developed a numerical method and correspond- 
ing computer code GALPROP for the calculation of Galactic 
CR propagation in 3D (Strong and Moskalenko, 1998). The 
code has been shown to reproduce simultaneously observa- 
tional data of many kinds related to CR origin and propaga- 
tion (Strong and Moskalenko, 1998; Moskalenko and Strong, 
1998; Moskalenko et al., 1998; Strong et al., 2000). The 
code has been validated on direct measurements of nuclei, 
antiprotons, electrons, and positrons, and astronomical mea- 
surements of 7-rays and synchrotron radiation. These data 
provide many independent constraints on model parameters. 

Here we use the GALPROP code for accurate calculation 
of production and propagation of secondary antiprotons. We 
explore the dependence of the antiproton flux on the nucleon 
injection spectrum and propagation parameters. The antipro- 
ton production is calculated using the pp production cross 
section and DTUNUC nuclear factors (Simon et al., 1998) 
or the pp production cross section scaled appropriately with 
atomic numbers. Inelastic scattering producing "tertiary" an- 
tiprotons and "secondary" protons is taken into account. The 
calculated local interstellar spectrum (LIS) is modulated us- 
ing the steady-state drift model. For the calculation reported 
here, we use a cylindrically symmetrical Galactic geometry. 



2 Propagation models and parameters 

The propagation parameters have been fixed using the B/C 
ratio. Nucleon injection spectra were chosen to reproduce 
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Table 1. Propagation parameter sets. 



Model 


Injection 
index", 7 


Diffusion coefficient'' 


Reacceleration/Convection 


Do, cm^ s ^ 


Index, 6 


WA.kms"^ dV/ dz, km kpc~^ 


Stochastic Reacceleration (SR) 


2.43 


6.10 X 10^* 


0.33 


30 


Minimal Reacceleration & Convection (MRC) 


2.43 


4.30 X 10^* 


0.33 


17 10 


Plain Diffusion (PD) 


2.16 


3.10 X 10^* 


0.6 




Diffusion plus Convection (DC) 


2.46/2.16'' 


2.50 X 10^® 


0./0.6'' 


10 



"For a power-law in rigidity, oc p 

= f3Do (p/4 GV)*, index S is shown below/above 4 GV. 
''Index below/above rigidity 20 GV. 



the local CR measurements. The source abundances of all 
isotopes Z < 28 are given in Strong and Moskalenko (2001). 
We thus use the same (Webber et al.) cross-section parametri- 
zation as in that work in our calculations. 

In all cases the halo size has been set to = 4 kpc, which 
is within the range = 3 — 7 kpc derived using the GAL- 
RPOP code and the combined measurements of radioactive 
isotope abundances, ^°Be, ^^Al, ^^Cl, and ^^Mn (Strong and 
Moskalenko 2001 and references therein). Note that the ex- 
act value of Zh is unimportant for antiproton calculations pro- 
vided that the propagation parameters are tuned to match the 
B/C ratio. 

The "tertiary" antiprotons (inelastically scattered secon- 
daries), significant at the lowest energies, are important in in- 
terstellar space, but make no difference when compared with 
measurements in the heliosphere. 

To investigate the range of interstellar spectra and prop- 
agation parameters we considered four basic models (Table 
1). Our results are plotted in Figs. 1-3. 

A model with the stochastic reacceleration (SR) reproduces 
the sharp peak in secondary to primary nuclei ratios in a 
physically understandable way without breaks in the diffu- 
sion coefficient and/or the injection spectrum (Strong and 
Moskalenko, 1998, 2001). However, this model produces a 
bump in proton and He spectra at ~ 2 GeV/nucleon which 
is not observed.' This bump can be removed by choosing an 
injection spectrum that hardens at low energies (Jones et al., 
2001). There are however some problems with secondaries 
such as positrons and antiprotons that are more difficult to 
manage. A similar bump appears in the positron spectrum at 
~ 1 GeV (Fig. 3), and the model underproduces antiprotons 
at 2 GeV by more than 30% (Fig. 2). Taken together they 
provide evidence against strong reacceleration^ in the ISM. 

Another model combining reduced reacceleration and con- 
vection (MRC) also produces too few antiprotons. 

Using a plain diffusion model (PD) we can get good agree- 
ment with B/C above few GeV/nucleon, with nucleon spectra 
and positrons, but this model overproduces antiprotons at 2 
GeV by ^ 20% and contradicts the secondary/primjiry nuclei 
ratio (B/C) below 1 GeV/nucleon. 



'a similar bump is produced also in the electron spectrum. 
^We define the reacceleration to be "strong" if the model is able 

to match the B/C ratio without invoking other mechanisms such as 
convection and/or breaks in the diffusion coefficient. 
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Fig. 1. Calculated proton LIS and modulated spectra ($ = 550 
MY). The line coding: solid lines - DC model (thick line - modu- 
lated spectrum), dashes - DR, dots - PD. The upper curve is always 
LIS spectrum, the lower is modulated. Data: IMAX (Menn et al., 
2000), CAPRICE (Boezio et al., 1999), AMS (Alcaraz et al., 2000), 
BESS (Sanuki et al., 2000). 



A model including diffusion and convection (DC) is our 
best fitting model. It reproduces all the particle data "on av- 
erage", although it has still some problem with the reproduc- 
tion of the sharp peak in B/C ratio. In this model a flattening 
of the diffusion coefficient below 4 GV is required to match 
the B/C ratio at low energies. 

To better match primaries (p, He) in the DC model, we 
introduced a steeper injection spectrum below 20 GV; such a 
break, however, has almost no effect on secondaries (p, e"*"). 
The existence of a sharp upturn below a few GeV/nucleon 
follows from SNR shock acceleration theory (Ellison et al., 
2001); this is a transition region between thermal and non- 
thermal particle populations in the shock. Our model does 
not require a sharp break (0.3 in index is enough). 

We use the DC model to calculate the LIS spectra of pro- 
tons and antiprotons and then use the drift model to deter- 
mine their modulated spectra and ratio over the solar cy- 
cles with positive (A > 0) and negative (A < 0) polarity 
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Fig. 2. Calculated spectrum of secondary p's for 3 models. For 
each model, the upper curve is the LIS, the lower curve - modulated 
(# = 550 MV). Line coding as in Fig. 1. The two lowest curves 
show separately the LIS and modulated spectrum of the "tertiary" 
component in DC model. Data: BESS (Orito et al., 2000), MASS91 
(Basini et al., 1999), CAPRICE98 (Boezio et al., 2001). 



(Figs. 4, 5). The variations shown depend on the tilt angle.-' 
This allows us to estimate the near-Earth spectra for arbitrary 
epochs in the past as well as make some predictions for the 
future. It may be also used to test the theory of heliospheric 
modulation (Langner and Potgieter, 2001). 



3 Discussion 

It appears quite difficult to get agreement with B/C, p & 
He spectra, and antiprotons spectrum simultaneously in the 
framework of simple "physical" reacceleration/convection 
models; this conclusion is mainly the result of the increased 
precision of the CR experimental data, and also the improved 
reliability of the calculations. 

What could be the origin of this failure, apart from the 
propagation models ? For B/C, it seems unlikely that the 
cross-sections are sigificantly in error since for B production 
from C and O they are well measured at the relevant ener- 
gies. For p and He the GeV bump produced by reaccelera- 
tion seems to be outside the limits allowed by the modulation 
models and the observed fluxes, although the uncertainties in 
modulation are still considerable. Such a bump in the spectra 
of primaries (p, He, e~) could be "removed" by including a 
flattening of the injection spectrum at low energies, while for 
secondaries (e+, p) we have no such freedom and the prob- 
lem remains. 

If we assume then that the problem is in the propagation 
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Fig. 3. Calculated positron flux for 3 models. For each model, the 
upper curve is the LIS, the lower curve - modulated ($ = 600 
MV). Line coding as in Fig. 1. Data: HEAT94 (Barwick et al., 
1998), CAPRICE94 (Boezio et al., 2000). 



models, we have shown that it is possible to construct a model 
(DC) which fits all these data, by postulating a significant 
flattening of the diffusion coefficient below 4 GV together 
with convection (and possibly with reduced reacceleration). 
This type of break in the diffusion coefficient is reminiscent 
of the standard procedure in "leaky-box" models where the 
escape time is set to a constant below a few GeV. This has 
always appeared a completely ad hoc device without phys- 
ical justification, but the present analysis suggest it may be 
forced on us by these type of data, so that possibiUties for its 
physical origin should be studied. At the same time the sen- 
sitivity of the whole analysis to the modulation models has 
to be investigated further. 

Acknowledgements. We thank Frank lones and Vladimir Ptuskin 
for fruitful discussions, and Don Ellison for discussions and a copy 
of his manuscript prior to publication. IVM acknowledges support 
from NAS/NRC Research Associateship Program. 



References 

Alcaraz, J., et al., Phys. Lett. B, 490, 27-35, 2000. 

Barwick, S. W., et al., Astrophys. I., 498, 779-789, 1998. 

Basini, G., et al., Proc. 26th ICRC (Salt Lake City), 3, 77-80, 1999. 

Bergstrom, D., et al., Astrophys. J., 534, L177-L180, 2000. 

Boezio, M., et al., Astrophys. J., 518, 457^72, 1999. 

Boezio, M., et al., Astrophys. I., 532, 653-669, 2000. 

Boezio, M., et al., astro-ph/0103513, 2001. 

Ellison, D. C, Slane, P., and Gaensler, B. M., Astrophys. J., sub- 
mitted, 2001. 

Jones, F. C, et al., Astrophys. J., 547, 264-271, 2001. 

Maeno, T, et al, astro-ph/0010381, 2001. 

Menn, W., et al., Astrophys. J., 533, 281-297, 2000. 



4 



T — I — I I I I I 1 1 1 — I — I I I I I 1 1 1 — I — I I I I I 1 — I I 1 — I — I I I I I 1 1 1 — I — I I I I I 1 1 1 — I — I I I I I r 




Kinetic energy, GeV Kinetic energy, GeV 

Fig. 4. Calculated proton LIS and modulated spectra for the two magnetic polarity dependent modulation epochs, A > (left) and A < 
(right). Tilt angle from top to bottom: 5° , 15° , 25° , 35° , 45° , 55° , 65° , 75° . The tilt angle corresponding to BESS and AMS data is ~ 5°- 15° 
{A > 0) depending on the coronal field model. On the right panel, ^ < 0, the data are shown only for guidance. Data refences as in Fig. 1. 
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Fig. 5. Calculated antiproton LIS and modulated spectra for the two magnetic polarity dependent modulation epochs, A > (left) and 
A<0 (right). Tilt angle from top to bottom: 5° , 15° , 25° , 35° , 45° , 55° , 65° , 75° . The tilt angle corresponding to BESS data is - 5° - 15° 
{A > 0) depending on the coronal field model. On the right panel, A < 0, the data are shown only for guidance. Data references as in Fig. 2. 
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